Validating saliva as a biological sample for cost-effective, rapid and routine screening for SARS-CoV-2.

PURPOSE: Compared to nasopharyngeal/oropharyngeal swabs (N/OPS-VTM), non-invasive saliva samples have enormous potential for scalability and routine population screening of SARS-CoV-2. In this study, we investigate the efficacy of saliva samples relative to N/OPS-VTM for use as a direct source for RT-PCR based SARS-CoV-2 detection. METHODS: We collected paired nasopharyngeal/oropharyngeal swabs and saliva samples from suspected positive SARS-CoV-2 patients and tested using RT-PCR. We used generalized linear models to investigate factors that explain result agreement. Further, we used simulations to evaluate the effectiveness of saliva-based screening in restricting the spread of infection in a large campus such as an educational institution. RESULTS: We observed a 75.4% agreement between saliva and N/OPS-VTM, that increased drastically to 83% in samples stored for less than three days. Such samples processed within two days of collection showed 74.5% test sensitivity. Our simulations suggest that a test with 75% sensitivity, but high daily capacity can be very effective in limiting the size of infection clusters in a workspace. Guided by these results, we successfully implemented a saliva-based screening in the Bangalore Life Sciences Cluster (BLiSC) campus. CONCLUSION: These results suggest that saliva may be a viable alternate source for SARS-CoV-2 surveillance if samples are processed immediately. Although saliva shows slightly lower sensitivity levels when compared to N/OPS-VTM, saliva collection is logistically advantageous. We strongly recommend the implementation of saliva-based screening strategies for large workplaces and in schools, as well as for population-level screening and routine surveillance as we learn to live with the SARS-CoV-2 virus.

Molecular and cellular constraints on vesicle traffic evolution.

In eukaryotic cells, the budding and fusion of intracellular transport vesicles is carefully orchestrated in space and time. Locally, a vesicle's source compartment, its cargo, and its destination compartment are controlled by dynamic multi-protein specificity modules. Globally, vesicle constituents must be recycled to ensure homeostasis of compartment compositions. The emergence of a novel vesicle pathway therefore requires new specificity modules as well as new recycling routes. Here, we review recent research on local (molecular) constraints on gene module duplication and global (cellular) constraints on intracellular recycling. By studying the evolution of vesicle traffic, we may discover general principles of how complex traits arise via multiple intermediate steps.

Genome doubling enabled the expansion of yeast vesicle traffic pathways.

Vesicle budding and fusion in eukaryotes depend on a suite of protein types, such as Arfs, Rabs, coats and SNAREs. Distinct paralogs of these proteins act at distinct intracellular locations, suggesting a link between gene duplication and the expansion of vesicle traffic pathways. Genome doubling, a common source of paralogous genes in fungi, provides an ideal setting in which to explore this link. Here we trace the fates of paralog doublets derived from the 100-Ma-old hybridization event that gave rise to the whole genome duplication clade of budding yeast. We find that paralog doublets involved in specific vesicle traffic functions and pathways are convergently retained across the entire clade. Vesicle coats and adaptors involved in secretory and early-endocytic pathways are retained as doublets, at rates several-fold higher than expected by chance. Proteins involved in later endocytic steps and intra-Golgi traffic, including the entire set of multi-subunit and coiled-coil tethers, have reverted to singletons. These patterns demonstrate that selection has acted to expand and diversify the yeast vesicle traffic apparatus, across species and time.

Graph-theoretic constraints on vesicle traffic networks.

Eukaryotic cells use small membrane-enclosed vesicles to transport molecular cargo between intracellular compartments. Interactions between molecules on vesicles and compartments determine the source and target compartment of each vesicle type. The set of compartment and vesicle types in a cell define the nodes and edges of a transport graph known as the vesicle traffic network. The transmembrane SNARE proteins that regulate vesicle fusion to target compartments travel in cycles through the transport graph, but the paths they follow must be tightly regulated to avoid aberrant vesicle fusion. Here we use graph-theoretic ideas to understand how such molecular constraints place constraints on the structure of the transport graph. We identify edge connectivity (the minimum number of edges that must be removed to disconnect a graph) as a key determinant that separates allowed and disallowed types of transport graphs. As we increase the flexibility of molecular regulation, the required edge connectivity decreases, so more types of vesicle transport graphs are allowed. These results can be used to aid the discovery of new modes of molecular regulation and new vesicle traffic pathways.

GRASP55 regulates intra-Golgi localization of glycosylation enzymes to control glycosphingolipid biosynthesis.

The Golgi apparatus, the main glycosylation station of the cell, consists of a stack of discontinuous cisternae. Glycosylation enzymes are usually concentrated in one or two specific cisternae along the cis-trans axis of the organelle. How such compartmentalized localization of enzymes is achieved and how it contributes to glycosylation are not clear. Here, we show that the Golgi matrix protein GRASP55 directs the compartmentalized localization of key enzymes involved in glycosphingolipid (GSL) biosynthesis. GRASP55 binds to these enzymes and prevents their entry into COPI-based retrograde transport vesicles, thus concentrating them in the trans-Golgi. In genome-edited cells lacking GRASP55, or in cells expressing mutant enzymes without GRASP55 binding sites, these enzymes relocate to the cis-Golgi, which affects glycosphingolipid biosynthesis by changing flux across metabolic branch points. These findings reveal a mechanism by which a matrix protein regulates polarized localization of glycosylation enzymes in the Golgi and controls competition in glycan biosynthesis.

Golgi compartments enable controlled biomolecular assembly using promiscuous enzymes.

The synthesis of eukaryotic glycans - branched sugar oligomers attached to cell-surface proteins and lipids - is organized like a factory assembly line. Specific enzymes within successive compartments of the Golgi apparatus determine where new monomer building blocks are linked to the growing oligomer. These enzymes act promiscuously and stochastically, causing microheterogeneity (molecule-to-molecule variability) in the final oligomer products. However, this variability is tightly controlled: a given eukaryotic protein type is typically associated with a narrow, specific glycan oligomer profile. Here, we use ideas from the mathematical theory of self-assembly to enumerate the enzymatic causes of oligomer variability and show how to eliminate each cause. We rigorously demonstrate that cells can specifically synthesize a larger repertoire of glycan oligomers by partitioning promiscuous enzymes across multiple Golgi compartments. This places limits on biomolecular assembly: glycan microheterogeneity becomes unavoidable when the number of compartments is limited, or enzymes are excessively promiscuous.

Promiscuity and specificity of eukaryotic glycosyltransferases.

Glycosyltransferases are a large family of enzymes responsible for covalently linking sugar monosaccharides to a variety of organic substrates. These enzymes drive the synthesis of complex oligosaccharides known as glycans, which play key roles in inter-cellular interactions across all the kingdoms of life; they also catalyze sugar attachment during the synthesis of small-molecule metabolites such as plant flavonoids. A given glycosyltransferase enzyme is typically responsible for attaching a specific donor monosaccharide, via a specific glycosidic linkage, to a specific moiety on the acceptor substrate. However these enzymes are often promiscuous, able catalyze linkages between a variety of donors and acceptors. In this review we discuss distinct classes of glycosyltransferase promiscuity, each illustrated by enzymatic examples from small-molecule or glycan synthesis. We highlight the physical causes of promiscuity, and its biochemical consequences. Structural studies of glycosyltransferases involved in glycan synthesis show that they make specific contacts with 'recognition motifs' that are much smaller than the full oligosaccharide substrate. There is a wide range in the sizes of glycosyltransferase recognition motifs: highly promiscuous enzymes recognize monosaccharide or disaccharide motifs across multiple oligosaccharides, while highly specific enzymes recognize large, complex motifs found on few oligosaccharides. In eukaryotes, the localization of glycosyltransferases within compartments of the Golgi apparatus may play a role in mitigating the glycan variability caused by enzyme promiscuity.

Frustration and fidelity in influenza genome assembly.

The genome of the influenza virus consists of eight distinct single-stranded RNA segments, each encoding proteins essential for the viral life cycle. When the virus infects a host cell, these segments must be replicated and packaged into new budding virions. The viral genome is assembled with remarkably high fidelity: experiments reveal that most virions contain precisely one copy of each of the eight RNA segments. Cell-biological studies suggest that genome assembly is mediated by specific reversible and irreversible interactions between the RNA segments and their associated proteins. However, the precise inter-segment interaction network remains unresolved. Here, we computationally predict that tree-like irreversible interaction networks guarantee high-fidelity genome assembly, while cyclic interaction networks lead to futile or frustrated off-pathway products. We test our prediction against multiple experimental datasets. We find that tree-like networks capture the nearest-neighbour statistics of RNA segments in packaged virions, as observed by electron tomography. Just eight tree-like networks (of a possible 262 144) optimally capture both the nearest-neighbour data and independently measured RNA-RNA binding and co-localization propensities. These eight do not include the previously proposed hub-and-spoke and linear networks. Rather, each predicted network combines hub-like and linear features, consistent with evolutionary models of interaction gain and loss.

How contraction has shaped evolution.

Two unicellular relatives of animals reveal that coordinated contractions of groups of cells using actomyosin predated animal multicellularity during evolution.

Correction: Population Distribution Analyses Reveal a Hierarchy of Molecular Players Underlying Parallel Endocytic Pathways.

[This corrects the article DOI: 10.1371/journal.pone.0100554.].

Small GTPases and BAR domain proteins regulate branched actin polymerisation for clathrin and dynamin-independent endocytosis.

Using real-time TIRF microscopy imaging, we identify sites of clathrin and dynamin-independent CLIC/GEEC (CG) endocytic vesicle formation. This allows spatio-temporal localisation of known molecules affecting CG endocytosis; GBF1 (a GEF for ARF1), ARF1 and CDC42 which appear sequentially over 60 s, preceding scission. In an RNAi screen for BAR domain proteins affecting CG endocytosis, IRSp53 and PICK1, known interactors of CDC42 and ARF1, respectively, were selected. Removal of IRSp53, a negative curvature sensing protein, abolishes CG endocytosis. Furthermore, the identification of ARP2/3 complex at CG endocytic sites, maintained in an inactive state reveals a function for PICK1, an ARP2/3 inhibitor. The spatio-temporal sequence of the arrival and disappearance of the molecules suggest a mechanism for a clathrin and dynamin-independent endocytic process. Coincident with the loss of PICK1 by GBF1-activated ARF1, CDC42 recruitment leads to the activation of IRSp53 and the ARP2/3 complex, resulting in a burst of F-actin polymerisation potentially powering scission.

Using stochastic cell division and death to probe minimal units of cellular replication.

The invariant cell initiation mass measured in bacterial growth experiments has been interpreted as a minimal unit of cellular replication. Here we argue that the existence of such minimal units induces a coupling between the rates of stochastic cell division and death. To probe this coupling we tracked live and dead cells in Escherichia coli populations treated with a ribosome-targeting antibiotic. We find that the growth exponent from macroscopic cell growth or decay measurements can be represented as the difference of microscopic first-order cell division and death rates. The boundary between cell growth and decay, at which the number of live cells remains constant over time, occurs at the minimal inhibitory concentration (MIC) of the antibiotic. This state appears macroscopically static but is microscopically dynamic: division and death rates exactly cancel at MIC but each is remarkably high, reaching 60% of the antibiotic-free division rate. A stochastic model of cells as collections of minimal replicating units we term 'widgets' reproduces both steady-state and transient features of our experiments. Sub-cellular fluctuations of widget numbers stochastically drive each new daughter cell to one of two alternate fates, division or death. First-order division or death rates emerge as eigenvalues of a stationary Markov process, and can be expressed in terms of the widget's molecular properties. High division and death rates at MIC arise due to low mean and high relative fluctuations of widget number. Isolating cells at the threshold of irreversible death might allow molecular characterization of this minimal replication unit.

Discovering vesicle traffic network constraints by model checking.

A eukaryotic cell contains multiple membrane-bound compartments. Transport vesicles move cargo between these compartments, just as trucks move cargo between warehouses. These processes are regulated by specific molecular interactions, as summarized in the Rothman-Schekman-Sudhof model of vesicle traffic. The whole structure can be represented as a transport graph: each organelle is a node, and each vesicle route is a directed edge. What constraints must such a graph satisfy, if it is to represent a biologically realizable vesicle traffic network? Graph connectedness is an informative feature: 2-connectedness is necessary and sufficient for mass balance, but stronger conditions are required to ensure correct molecular specificity. Here we use Boolean satisfiability (SAT) and model checking as a framework to discover and verify graph constraints. The poor scalability of SAT model checkers often prevents their broad application. By exploiting the special structure of the problem, we scale our model checker to vesicle traffic systems with reasonably large numbers of molecules and compartments. This allows us to test a range of hypotheses about graph connectivity, which can later be proved in full generality by other methods.

Organelle acidification: an ancient cellular leak detector.

Intracellular membrane-bounded organelles of eukaryotic cells transiently contact the extracellular environment during endocytosis and secretion. Such contacts must be precisely timed to prevent leakage of cargo. I argue that early eukaryotes evolved organelle acidification as a way to detect and prevent leakage.

Stacking the odds for Golgi cisternal maturation.

What is the minimal set of cell-biological ingredients needed to generate a Golgi apparatus? The compositions of eukaryotic organelles arise through a process of molecular exchange via vesicle traffic. Here we statistically sample tens of thousands of homeostatic vesicle traffic networks generated by realistic molecular rules governing vesicle budding and fusion. Remarkably, the plurality of these networks contain chains of compartments that undergo creation, compositional maturation, and dissipation, coupled by molecular recycling along retrograde vesicles. This motif precisely matches the cisternal maturation model of the Golgi, which was developed to explain many observed aspects of the eukaryotic secretory pathway. In our analysis cisternal maturation is a robust consequence of vesicle traffic homeostasis, independent of the underlying details of molecular interactions or spatial stacking. This architecture may have been exapted rather than selected for its role in the secretion of large cargo.

On the Archaeal Origins of Eukaryotes and the Challenges of Inferring Phenotype from Genotype.

If eukaryotes arose through a merger between archaea and bacteria, what did the first true eukaryotic cell look like? A major step toward an answer came with the discovery of Lokiarchaeum, an archaeon whose genome encodes small GTPases related to those used by eukaryotes to regulate membrane traffic. Although 'Loki' cells have yet to be seen, their existence has prompted the suggestion that the archaeal ancestor of eukaryotes engulfed the future mitochondrion by phagocytosis. We propose instead that the archaeal ancestor was a relatively simple cell, and that eukaryotic cellular organization arose as the result of a gradual transfer of bacterial genes and membranes driven by an ever-closer symbiotic partnership between a bacterium and an archaeon.

Wine glasses and hourglasses: Non-adaptive complexity of vesicle traffic in microbial eukaryotes.

Microbial eukaryotes present a stunning diversity of endomembrane organization. From specialized secretory organelles such as the rhoptries and micronemes of apicomplexans, to peroxisome-derived metabolic compartments such as the glycosomes of kinetoplastids, different microbial taxa have explored different solutions to the compartmentalization and processing of cargo. The basic secretory and endocytic system, comprising the ER, Golgi, endosomes, and plasma membrane, as well as diverse taxon-specific specialized endomembrane organelles, are coupled by a complex network of cargo transport via vesicle traffic. It is tempting to connect form to function, ascribing biochemical roles to each compartment and vesicle of such a system. Here we argue that traffic systems of high complexity could arise through non-adaptive mechanisms via purely physical constraints, and subsequently be exapted for various taxon-specific functions. Our argument is based on a Boolean mathematical model of vesicle traffic: we specify rules of how compartments exchange vesicles; these rules then generate hypothetical cells with different types of endomembrane organization. Though one could imagine a large number of hypothetical vesicle traffic systems, very few of these are consistent with molecular interactions. Such molecular constraints are the bottleneck of a metaphorical hourglass, and the rules that make it through the bottleneck are expected to generate cells with many special properties. Sampling at random from among such rules represents an evolutionary null hypothesis: any properties of the resulting cells must be non-adaptive. We show by example that vesicle traffic systems generated in this random manner are reminiscent of the complex trafficking apparatus of real cells.

Universal Poisson Statistics of mRNAs with Complex Decay Pathways.

Messenger RNA (mRNA) dynamics in single cells are often modeled as a memoryless birth-death process with a constant probability per unit time that an mRNA molecule is synthesized or degraded. This predicts a Poisson steady-state distribution of mRNA number, in close agreement with experiments. This is surprising, since mRNA decay is known to be a complex process. The paradox is resolved by realizing that the Poisson steady state generalizes to arbitrary mRNA lifetime distributions. A mapping between mRNA dynamics and queueing theory highlights an identifiability problem: a measured Poisson steady state is consistent with a large variety of microscopic models. Here, I provide a rigorous and intuitive explanation for the universality of the Poisson steady state. I show that the mRNA birth-death process and its complex decay variants all take the form of the familiar Poisson law of rare events, under a nonlinear rescaling of time. As a corollary, not only steady-states but also transients are Poisson distributed. Deviations from the Poisson form occur only under two conditions, promoter fluctuations leading to transcriptional bursts or nonindependent degradation of mRNA molecules. These results place severe limits on the power of single-cell experiments to probe microscopic mechanisms, and they highlight the need for single-molecule measurements.

Ancient dynamin segments capture early stages of host-mitochondrial integration.

Eukaryotic cells use dynamins-mechano-chemical GTPases--to drive the division of endosymbiotic organelles. Here we probe early steps of mitochondrial and chloroplast endosymbiosis by tracing the evolution of dynamins. We develop a parsimony-based phylogenetic method for protein sequence reconstruction, with deep time resolution. Using this, we demonstrate that dynamins diversify through the punctuated transformation of sequence segments on the scale of secondary-structural elements. We find examples of segments that have remained essentially unchanged from the 1.8-billion-y-old last eukaryotic common ancestor to the present day. Stitching these together, we reconstruct three ancestral dynamins: The first is nearly identical to the ubiquitous mitochondrial division dynamins of extant eukaryotes, the second is partially preserved in the myxovirus-resistance--like dynamins of metazoans, and the third gives rise to the cytokinetic dynamins of amoebozoans and plants and to chloroplast division dynamins. The reconstructed sequences, combined with evolutionary models and published functional data, suggest that the ancestral mitochondrial division dynamin also mediated vesicle scission. This bifunctional protein duplicated into specialized mitochondrial and vesicle variants at least three independent times--in alveolates, green algae, and the ancestor of fungi and metazoans-accompanied by the loss of the ancient prokaryotic mitochondrial division protein FtsZ. Remarkably, many extant species that retain FtsZ also retain the predicted ancestral bifunctional dynamin. The mitochondrial division apparatus of such organisms, including amoebozoans, red algae, and stramenopiles, seems preserved in a near-primordial form.